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The decay channel J/ip -^ ^ujuj, uj -^ Tr+vr^vr'^ is analyzed using a sample of 5.8 
X 10 "^ J /ij) events collected with the BESII detector. The ujuj invariant mass dis- 
tribution peaks at 1.76 GeV/c^, just above the ujuj threshold. Analysis of angular 



correlations indicate that the ujoj system below 2 GeV/c^ is predominantly pseu- 
doscalar. A partial wave analysis confirms the predominant pseudoscalar structure, 
together with small 0^+ and 2++ contributions, and yields a pseudoscalar mass M = 
1744 ± 10 (stat) ± 15 (syst) MeV/c^, a width T = 244+^4 (gtat) ± 25 (syst) MeV/c^, 
and a product branching fraction Br(J/V' -^ 7ry(1760)) • Br(ry(1760) -^ lolo) = (1.98 
± 0.08 (stat) ± 0.32 (syst)) x 10^^. 

PACS numbers: 12.39. Mk, 13.20.Gd, 13.30.Ce, 14.40. Cs 

I. INTRODUCTION 

Glueballs are expected to be copiously produced in radiative J /ip decays. How- 
ever, until now, no unique experimental signatures of such states have been found. 
The pseudoscalar ground state mesons (1^5*0) are well established, and 7r(1300), 
ry(1295), r7(1475), and 7^(1460) are suggested as the first radial excitations (2^Sq) 
of the pseudoscalar mesons [1]. In the Particle Data Group (PDG) listings, two 
pseudoscalar states are reported in the r7(1440) mass region. However, there are 
too many pseudoscalar states, and it is very difficult to find a place for the lower 
mass 7^(1440) or the ?7(1760) within any qq model [1]. At one time, the ?7(1440) 
was regarded as a glueball candidate when it was observed in J /ip radiative de- 
cay [2] and there was only an upper limit on its two-photon production [3]. But 
this viewpoint changed when its radiative decay modes [4-7] were observed and it 
was also observed in untagged 77 collisions by the L3 collaboration [8] . In addition, 
lattice gauge theory would have great difficulty to accommodate such a low-mass 
0-+ glueball [9]. 

The 77(1 760) was reported by the MARK HI collaboration in J /ip radiative decays 
and was found to decay to uouo [10] and pp [11]. It was also observed by the DM2 
collaboration in J /%l) radiative decays in the pp decay mode with a mass of M = 
1760 ± 11 MeV/c^ and a width of T = 60 ± 16 MeV/c^ [12] and in the ujuj decay 
mode [13]. The BESI experiment reported its rjii^Ti^ decay with a mass of M = 
1760 ± 35 MeV/c^, but without a determination of its width [14]. Also, possible 
pseudoscalar production at threshold in the (fxp mode has been observed in tt~p 



scattering [15]. The ri{1760) was suggested to be a 3^5*0 pseudoscalar qq meson, but 
some authors suggest a mixture of glueball and qq or a hybrid [16, 17]. Recently, 
in Ref. [18], it was argued that the pseudoscalar glueball may be in the 1.5 to 1.9 
GeV/c^ mass region, and that it also has Vector Vector decay modes. In this paper, 
we present results from an analysis of J/ip -^ '-/ujuj , uj -^ ti^ti^ti^ decays, based on 
a sample of 58 million J/ip events collected with the BESII detector at the Beijing 
Electron- Positron Collider (BEPC). The presence of a signal around 1.76 GeV/c^ 
and its pseudoscalar character are confirmed, and the mass, width, and branching 
fraction are measured by partial wave analysis. 

II. BES DETECTOR AND MONTE CARLO SIMULATION 

BESII is a large solid-angle magnetic spectrometer that is described in de- 
tail in Ref. [19]. Charged particle momenta are determined with a resolution of 
(Tp/p = 1.78%vT+~p^ (with p in GeV/c) in a 40-layer cylindrical main drift cham- 
ber (MDC). Particle identification is accomplished using specific ionization (dE/dx) 
measurements in the MDC and time-of-fiight (TOP) measurements in a barrel-like 
array of 48 scintillation counters. The dE/dx resolution is CdE/dx = 8.0%; the TOP 
resolution is (Ttof = 180 ps for the Bhabha events. Outside of the TOP counters is a 
12-radiation-length barrel shower counter (BSC) comprised of gas tubes interleaved 
with lead sheets. The BSC measures the energies and directions of photons with 



resolutions of as/E ~ 21%/JE(GeV), a^ = 7.9 mrad, and az = 2.3 cm. The iron 
flux return of the magnet is instrumented with three double layers of counters that 
are used to identify muons. 

In this analysis, a GEANT3 based Monte Carlo (MC) simulation program 
(SIMBES) [20] with detailed consideration of real detector responses (such as dead 
electronic channels) is used. The consistency between data and Monte Carlo has 
been carefully checked in many high-purity physics channels, and the agreement is 
quite reasonable [20]. 



III. EVENT SELECTION 

J /ip ^ 7 + 2(7r+7r^7r°) candidates are selected from events with four charged 
tracks in the drift chamber and five photons in the barrel shower counter. 

A. Charged particle identification 

Each charged track, reconstructed using MDC information, is required to be 
well fitted to a helix, to be within the polar angle region | cos^| < 0.85, to have a 
transverse momentum larger than 50 MeV/c, and have the point of closest approach 
of the track to the beam axis within 2 cm of the beam axis and within 20 cm from 
the center of the interaction region along the bean line. For each track, we make 
a weak particle identification requirement: either the TOF or dE/dx information 
must agree with that expected for a pion within four standard deviations. 

B. Photon identification 

Each candidate photon is required to have an energy deposit in the BSC greater 
than 35 MeV, to be isolated from charged tracks by more than 6°, to have the angle 
between the cluster development direction in the BSC and photon emission direction 
less than 30°, and to have the first hit in the beginning six radiation lengths. 

C. Event selection criteria 

Events are required to have four charged tracks with net charge zero and 
have from 5 to 8 photon candidates. Six-constraint (6-C) kinematic fits to the 
J /ip — > 7 + 2(7r+7r^7r°) hypothesis are made with both 77 invariant masses be- 
ing constrained to the 7r° mass using all possible photon combinations. Note that 
there are 15 possible ways of combining five photons to obtain two 7r°s. We select the 
combination with the highest probability and require this probability to be greater 
than 10%. Six-constraint kinematic fits also are applied using the J /ip -^ 2(7r+7r^7r°) 
hypothesis, and the probability of these fits is required to be less than that of the 
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signal hypothesis. Two further requirements are imposed on events with more than 
five gammas to reduce the background from J/ip -^ 7r°2(7r+7r^7r'^). First, all gam- 
mas which do not belong to the chosen combination must have E^ < 140 MeV. 
Second, seven-constraint kinematic fits are performed to the J/ip -^ 7r'^2(7r^7r^7r°) 
hypothesis with the invariant mass of the three 77 pairs being constrained to the 
7r° mass, and the event is discarded if P(x^)7c > -P(x^)6c- 

The 7r"'"7r~7r° invariant mass distribution for the selected events is shown in Fig. 
1(a), where there are 8 entries per event. A clear u signal is present, mainly due to 
J/ip -^ WK^n^'n-^TT^. The open histogram in Fig. 1(b) shows the tt'^tt'tc^ invariant 
mass spectrum after the requirement that the invariant mass of the other tt^ti^tc^ 
is inside the u region (|m^+7r-7r0 — "^^1 < 40 MeV/c^), and the shaded histogram is 
for events after the requirement that the invariant mass of the other tt+tt^tt^ mass 
is inside the uj sideband region (40 MeV/c^ < |m^+^-7r0 — m^\ < 80 MeV/c^). 
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FIG. 1: (a) The tt^tt^tt^ invariant mass distribution (8 entries per event), (b) The vr'^vr^Tr'' 
mass distribution for the best coco combination after the requirement that the invariant 
mass of the other tt^tt^tt^ is inside the to region (open histogram), defined by \m^+^-^o — 
m^\ < 40 MeV/c^, or in the sideband range (shaded histogram), defined by 40 MeV/c^ < 
\m^+^-^o — m^\ < 80 MeV/c^. (c) The vr+vr^vr'^ versus the vr^vr^vr'^ invariant mass (4 
entries per event). 

The tt+tt^tt'^ versus tt+tt^tt^ invariant mass distribution (four entries per event) is 
plotted in Fig. 1(c). A cluster of events is observed corresponding to uu production. 
Because the processes J/tp —^ ujuj and J /ip -^ -k^ujuj are forbidden by C-invariance, 
the presence of two tu's is direct evidence for the radiative decay J /ip — > ''yuju. The 



histogram of Fig. 2(a) shows the 2(7r+7r^7r°) invariant mass distribution of events 
with both tt+tt^tt'^ masses within the uj range (|m7r+7r-7rO —'m^\ <40 MeV/c^). There 
are 3046 events with a clear peak at 1.76 GeV/c^. The phase space invariant mass 
distribution and the acceptance versus ujuj invariant mass are also shown in the 
figure. The corresponding Dalitz plot is shown in Fig. 2(b). 
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FIG. 2: (a) The 2(7r+7r^7r'^) invariant mass distribution for candidate events. The dashed 
curve is the phase space invariant mass distribution, and the dotted curve shows the 
acceptance versus the cjo; invariant mass, (b) The Dalitz plot, (c) The 2(7r"*'7r~7r'^) invariant 
mass of the inclusive Monte Carlo sample (shaded histogram) . 



From MC simulation, the u signal can be well fitted with a double Gaussian 
with widths 14.6 MeV/c^ and 45.5 MeV/c^, which makes it difficult to evaluate the 
background from the sideband events especially for the events near ujuj threshold. 
To ensure that the structure at the ujuj mass threshold is not due to background, we 
have made studies of potential background processes. The main background sources 
come from J/ip -^ ujti^ti^ti^tt^, 72(7r+7r^7r°), and 7r°2(7r+7r^7r°). More than one half 
of the background comes from the first case. However, none of these background 
channels gives a peak at 1.76 GeV/c^ in the 2(7r+7r^7r°) invariant mass spectrum. 
In addition, possible backgrounds were checked with a MC sample of 58 million 
inclusive J ftp decays generated by the LUND model [21], with 'jujuj events removed. 
The shaded histogram of Fig. 2(c) shows the 2(7r+7r^7r'^) invariant mass distribution 
of the inclusive sample. There is no peak at the ujuj mass threshold in the invari- 
ant mass distribution at around 1.76 GeV/c^, while the inclusive MC distribution is 



comparable with data for 2(7r+7r^7r°) invariant mass greater than 2 GeV/c^. A back- 
ground evaluation is performed by fitting the tt+tt^tt^ mass distribution for events 
with the other tt+tt^tt'' within the uo signal region (Fig. 1(b)). The background 
shape is obtained from the inclusive MC sample, and the signal shape is obtained 
from the phase space MC sample of J /ip — > '-yuouo. The number of events is free in 
the fitting. The fitting yields 1441 ± 50 background events within the ujuj invariant 
mass range from 1.6 GeV/c^ to 2.8 GeV/c^. 

IV. ANGULAR CORRELATION ANALYSIS 

An analysis of the angular distributions of the accepted events has been per- 
formed in order to estimate whether the ujoo production below 2 GeV/c^ belongs 
to a resonant state with definite spin-parity. Candidate events and side-band back- 
ground events are analyzed choosing the ujuj pair whose quadratic sum of the two 
differences (rriTr+Tr-TrO ~f^u}) is minimum. For systems of two vector mesons, the dis- 
tribution of Xi the azimuthal angle between the normals to the two uj decay planes 
in the ujuj rest frame, provides a unique signature for the spin and parity [22-25]. 
The distribution takes the form dN/dx oc 1 +/9 cos(2x), where /3 is a constant which 
is independent of the polarization of the ujuj system, but exhibits strong correlation 
with the spin-parity. (3 is zero for odd spin and non-zero for even spin. Its sign is 
the parity of the uuj system. For 7-^ = 0", where {3 is —1, dN/dx oc sin^ Xi ^"^^ the 
effect is maximal. 

Fig. 3 shows the x distribution for a) J /ip — > ''jujuj events with m^^ less than 
2 GeV/c^, and b) events from the uj sidebands, where uj sideband is defined with 
a least one mT,+T,-T,o mass in the range of 40 MeV/c^ < \mT,+T,-T,o — m^\ < 120 
MeV/c^. The distributions for signal and background events are strikingly different, 
indicating a large component with even spin and odd parity in the signal region. 
The solid line in Fig. 3(a) is the result of a fit to a + 6sin^ x, which yields a sin^ x 
contribution of the '-jujuj event candidates below 2 GeV/c^ of 38.3 ± 3.5%. 
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FIG. 3: Distribution of X) the azimuthal angle between the normals to the two lo decay 
planes for (a) for uj — uj signal events {m^^^ < 2 GeV/c^) and (b) for lo sideband events 
(m^^ < 2 GeV/c^) 

V. PARTIAL WAVE ANALYSIS 



A partial wave analysis (PWA) has been carried out for events with 2(7r+7r^7r°) 
invariant mass from 1.6 GeV/c^ to 2.8 GeV/c^. The sequential decay process can 
be described by J/tp — > '~fX,X — > ujuj and uj —>■ tt+tt^tt*^. The amplitudes of the 
two body or three body decays are constructed using the covariant helicity coupling 
amplitude method [26]. The intermediate resonance X is denoted with the normal 
Breit-Wigner propagator BW = 1/(M^ — s — IMF), where s is the ujuj invariant 
mass-squared and M and F are the resonance's mass and width. The amplitude of 
sequential decay process is the product of all decay amplitudes and the Breit-Wigner 
propagator. 

The X angular distribution shows a strong contribution from structures with even 
spin and odd parity for ujuj invariant mass below 2 GeV/c^. Therefore the study of 
the 77(1 760) is the main goal of this analysis. From the PDG, three /2 resonances 
can decay into ujuj final states, /2(1560), /2(1640), and /2(1910). Because the mass 
of /2(1560) and /2(1640) are very close, only one resonance, /2(1640), is considered 
in the analysis. The /o(1710) is a well known resonance, and its spin-parity allows 
it to decay to a ujuj final state, so it is included in the fit. Finally, four possible 
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intermediate resonances 77(1760), /o(1710), /2(1640), and /2(1910) are included in 
the final analysis, and the total differential cross section da/d^ is 

^ = \A{r^) + A{fo) + A{f^) + A(/|)p + BG, (1) 

where ^(77) , A(/o), ^(/2), and A(/|) are the total amplitudes of the resonances 
77(1760), /o(1710), /2(1640), and /2(1910), respectively, and BG denotes the back- 
ground contribution, which is described by phase space. 

The relative magnitudes and phases of the amplitudes are determined by an 
unbinned maximum likelihood fit. The basis of likelihood fitting is calculating the 
probability that a hypothesized probability distribution function would produce the 
data set under consideration. The joint probability density for observing the N 
events in the data sample is 

N 

C = l[Pix,) (2) 

j=i 

where P{xi) is the probability to produce event i characterized by the measurement 

Xi, which is the normalized differential cross section: 

Pf^^) = A^ (3) 



The normalization integral / -^d^ is done by the phase space MC sample, the 
details are described in Ref. [27]. The free parameters are optimized by MINUIT 
[28]. Technically, rather than maximizing C, the S = —InC is minimized. In the 
minimization procedure, a change in log likelihood of 0.5 represents a one standard 
deviation effect for one parameter case. 

For the production of a pseudoscalar, only V waves are allowed in both the 
radiative decay J/ip -^ 7X and the hadronic decay X -^ uuo. For the production of 
a scalar, both S and V waves are possible in both the radiative and hadronic decays, 
but only S wave is considered in the fit. For the production of a 2^ resonance, S 
waves in both decays are considered, and two of three V waves in the radiative decay 
and only one T> wave in the hadronic decay, corresponding to the lower overall spin 
of the uouo system, are considered. From the analysis of angular correlations, it is 
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found that the contributions from /o(1710), /2(1640), and /2(1910) are very small, so 
the mass and width of these resonances are fixed to PDG values, but the amplitudes 
are allowed to vary in the fit. The mass and width of the ri{1760) are obtained from 
the optimization; the mass and width are M = 1744 ±10 MeV/c^ and T = 244^21 
MeV/c^, where the errors are statistical. The final global fit, and the contributions 
of all resonances and backgrounds are shown in Fig. 4. 
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FIG. 4: The 2(7r"'"7r~7r^) invariant mass distribution for J/tjj -^ 'jlouj. The points with error 
bars are data, the full histograms show the projection of the maximum likelihood fit, and 
the dashed histograms show the contributions of each of the resonances and background. 



Comparisons of angles of fit projections and data are shown in Fig. 5. To 
determine the goodness of fit, a x^ is calculated by comparing the data and fit 
projection histograms, where x^ is defined as [27]: 
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FIG. 5: Comparisons of angular distributions between data and fit projections of the 
global fit. The dashed histograms show the background contributions, (a, b) The polar 
and azimuthal angles of the radiative gamma, (c, d) The polar and azimuthal angles of 
the to in the ujuj rest system, (e) The polar angle of the normal to lo decay plane in the uj 
rest system, both oj's angle are projected in the same plot, (f) The x distribution. 



and rii and Vi are the number of events for data and the fit projections in the ith 
bin of each figure, respectively. The x^ ^''^^ number of degrees of freedom (ndf) for 
the UJUJ invariant mass and the angle distributions are shown in Table I, where the 
number of bins is taken as the number of degrees of freedom. The values of 'x^ /ndf 
are between 0.6 and 1.8, indicating good agreement between data and the fit. 

The numbers of events, detection efficiencies, and the corresponding branching 
fractions for J/ip -^ 7X — * '-jujuj with intermediate resonances r7(1760), /o(1710), 
/2(1640), and /2(1910) are shown in the Table II, where the errors are statistical 
errors only and the correlations between the different resonances are included. The 
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TABLE I: Goodness of fit check for the invariant mass distribution and angular distribu- 
tions shown in Fig. 5, where ndf and C.L. are the number of degrees of freedom and the 
corresponding confidence level. 

magnitudes and phases of the partial amplitudes from the PWA are used in the de- 
tection efficiency determination. Details of the fitting procedure and the detection 
efficiency determination can be found in Ref. [27]. The changes of the log likeli- 
hood value AS when the corresponding resonance is dropped from the fit and the 
statistical significance for each component are also shown in Table II, where the 
significance is calculated from the difference between S values of the fits with and 
without the resonance. The product branching fraction is Br(J/-?/' -^ 777(1760)) ■ 
Br(r7(1760) -^ ujlo) = (1.98 ± 0.08 (stat)) x 10"^, and the statistical significance 
of the ?7(1760) is above 10 a. All the resonances listed improve the fitting by more 
than 5a. If the spin-parity of ?7(1760) is replaced by 0++ in the fit, the log likelihood 
is worse by 248.0, so the possibility that its spin-parity is O"*""*" is excluded by at least 
10 a. 

The fit determines 1371 ± 45 background events, which is consistent with the 
result obtained from Fig. 1(b). Another technique for treating background, which 
was used in Refs. [27, 29, 30], is to set BG to in Eq. (1) and to cancel the 
background contribution by including MC data in the fit with the opposite sign of 
log likelihood compared to the data. As a check we have also used this method. The 
MC sample is obtained using the inclusive MC, and the number of background events 
is fixed to the fitting results of Fig. 1(b). The ri{1760) mass and width obtained are 
1742 ± 10 MeV/c^ and 234 ± 17 MeV/c^, respectively, and the product branching 
fraction of J/^ ^-fX,X ^ uu is (2.01 ± 0.08) x 10~l 
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resonance 


Events 


eff ( % ) 


Sr(xl0-3) 


Sys Err( % ) 


AS 


Sig. 


r?(1760) 


1045 ± 41 


1.15 


1.98 ± 0.08 


16.4 


280 


> 10a 


/o(1710) 


180 ± 37 


1.27 


0.31 ±0.06 


25.1 


23.5 


6.5(7 


/2(1910) 


151 ±32 


1.68 


0.20 ± 0.04 


64.9 


23.5 


5.8(7 


/2(1640) 


141 ± 26 


1.08 


0.28 ±0.05 


59.6 


21.4 


5.5(7 



TABLE II: The fitted number of events, detection efficiency, product branching fraction, 
systematic error, log likelihood value differences, and statistical significance of each reso- 



nance. 



Other states listed in the PDG between 1.6 GeV/c^ and 2.0 GeV/c^, that are con- 
sistent with decay into uju under spin-parity constraints, are the 772 (1645), ri{1870), 
/2(1810), /2(1950). If these resonances are included in the fit, the log likelihood 
values S improve by 6.2, 9.0, 8.9, and 8.4, respectively, while the ?7(1760) masses, 
widths, and branching fractions are consistent with the final fit result within statis- 
tical errors. The difference between results including and not including the ri{1870) 
will be taken as a systematic error. 

The /o(1790) has been recently claimed in J/ip decay [29]. If the parameters of 
/o(1710) are replaced with those of the /o(1790), the log likelihood value is improved 
by 2.8 after the reoptimization. The reoptimized mass, width, and product branch- 
ing ratio of ?7(1760) are 1744±10 MeV/c^, 238±20MeV/c^ and (1.97±0.07) x lO'^, 
respectively, and the product branching ratio of /o is then (0.39 ±0.07) x 10^^ (sta- 
tistical error only). Recently a scalar enhancement, the /o(1812), near uxf) threshold 
in J/ip -^ '-jujcj) decay was reported by the BESII collaboration [31], and information 
on the corresponding oou decay mode is very important to understand its nature [32- 
37]. If the /o(1710) parameters are replaced with the /o(1812) [31], the log likelihood 
value is improved by 0.8 after the reoptimization. The reoptimized mass, width, and 
product branching ratio of ?7(1760) are 1740 ± 10 MeV/c^, 246 ± 24 MeV/c^, and 
(1.97 ± 0.07) X 10^^, respectively, and the product branching ratio of /o is then 
(0.26 ± 0.05) X lO"''^ (statistical error only). If both the /o(1710) and /o(1812) are 
added in the fitting, the log likelihood values is improved by 5.2. If the /o(1710)'s 
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parameters are replaced with those of the /o(2020), the log likelihood value will be 
improved by 2.9. But if no scalar in this energy region is used in the fit, the log 
likelihood value is worse by 23.5, corresponding to 6.5 a. From these tests, we con- 
clude that a scalar is needed, but it is very difficult to determine its mass and width 
accurately due to the dominant contribution of the pseudoscalar. If the parameters 
of the /2(1640) are replaced with those of /2(1560), the likelihood value is improved 
by 0.55 after mass and width reoptimization. If the parameters of the /2(1910) are 
replaced by those of the /2(1950), the log likelihood value is improved by 2.1. In the 
final fit, the uj decay amplitude is described with sequential two body decays with 
V wave. If the u decay amplitude is taken to be constant, the results do not change 
much. In all these tests, the 77(1760) masses, widths, and the branching fractions are 
consistent with the final fit results. The differences are included in the systematic 
errors. 

VI. SYSTEMATIC ERROR 

The systematic errors are estimated by considering the following: (a) The 
/o(1710) is replaced with the /o(2020). (b) The /2(1640) is replaced with the 
/2(1560). (c) The /2(1910) is replaced with the /2(1950). (d) The fit is done with 
and without the 772 (1870). (e) A constant u amplitude is used in the fit. (f) Differ- 
ent background treatments, (g) Different 7 selection criteria: energy greater than 
50 MeV/c^, and the minimum angle between the gamma and the nearest charged 
track greater than 10°. (h) Changing the polar angle requirement of charged tracks 
to |cos^| < 0.8. (i) Changing the 6-C kinematic fit probability requirement from 
Probgc > 0.1 to Probgc > 0.05. (j) Changing the tt^tt^tt^ invariant mass requirement 
from the 40 MeV/c^ to 45 MeV/c^. The total errors are obtained by adding the indi- 
vidual errors in quadrature. The total mass and width systematic errors are 0.83% 
and 10.5%, respectively. For the branching fraction systematic error, the uncer- 
tainties in the MDC tracking, the photon identification efficiency, the 00 -^ tt+tt^tt^ 
branching fraction, and the number of J/ip events are also included, and the total 
branching fraction systematic errors are 16.4%, 25.1%, 64.9%, and 59.6% for 7^(1760), 
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/o(1710), /2(1910), and /2(1640), respectively, which are also listed in Table II. 

VII. DISCUSSION 

The ?7(1760) is prominently produced in J /ip — > '-yujuo. Its two-gluon coupling can 
discriminate between its gluonic and qq nature [16, 17]. If perturbative QCD works 
well and the non-relativistic approximation is applicable, the formalism proposed 
in Refs. [38, 39], which connects the two-gluon width r(r7(1760) -^ gg) of ?7(1760) 
to the radiative J/ip branching fraction, Br{J/il) -^ 77^(1760)), can be used (see 
Ref. [38], Eq. 3.4): 

where M is the r7(1760) mass, rj,(i76o)^3(, i^ ^^^ width to gg, and x\Hps{x)\'^ is 
the magnitude of the loop integral calculated in Ref. [38]. Therefore, the glu- 
onic content of ?7(1760) is estimated by its two-gluon coupling, which is calculated 
from its mass, width, and branching fraction in J/ip radiative decay. Rewriting 
^i]{i760)-^gg ^s rj,(i76o) ■-8^(^(1760) -^ gg) and multiplying both sides of the equation 
by Br{ri{1760) -^ ujuj), we obtain: 



fir(?7(1760) -^ gg) ■ fir(?7(1760) -^ cucu) 
= 10^ [BriJ/ij -^ 7?7(1760)) ■ Briri{l760) -^ ucu)] 



sr^/w, ,„,,„,,„, ,, , 1.5 GeV/c2\ /50 MeV/c2\ 45 

M l\ r^(i760) ; 39 



98 ± 16MeV _ o.os 

244tiMeV " ^'^^-o-o^' 



where x\Hps{x)\'^ is taken as 39, which is obtained from Fig.l of Ref. [38], and the 
theoretical uncertainty is not considered. Since we expect i?r(?7( 1760) -^ uouo) < 1.0, 
the relationship above implies -Br(r7(1760) -^ gg) > 0.28 at 90% confidence level. 

From a theoretical viewpoint, the coupling of a glueball to two photons is expected 
to be very weak, so the study of 77(1 760) production in the two photon process 
is needed. The 7^(1760) is abundantly produced in the J/ip radiative decay, but 
it is not seen in J/ip -^ '~f^V{p,(f)) [4-7], which means that the partial width of 
7^(1760) -^ 'yV{p, 0) is very small. ?7(1760) is shown to have large gluon component. 
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but its mass is much lower than the prediction from lattice QCD calculation [9], 
suggesting that it is a mixture of the glueball and qq meson. If r7(1760) is a mixed 
pseudoscalar glueball candidate, it should have flavor symmetric decays. Therefore 
other decay modes of r7(1760), such as ?7(1760) -^ pp, K*K*, rjinT, KKtt, etc. should 
be studied. Further studies are needed to understand the nature of the ri{1760), both 
experimentally and theoretically. 

VIII. SUMMARY 

In summary, J/ip -^ '-/ujuj, oo -^ 7r+7r^7r° is studied, and the ujoo invariant mass dis- 
tribution peaks at 1.76 GeV/c^. The partial wave analysis shows that the structure 
is predominantly pseudoscalar, with small contributions from /o(1710), /2(1640), 
and /2(1910). The mass of the pseudoscalar is M = 1744 ± 10 (stat) ± 15 (syst) 
MeV/c^ the width T = 244lii (stat) ± 25 (syst) MeV/c^ and the product branch- 
ing fraction is Br(J/^ -^ 77^(1760)) ■ Br(r7(1760) -^ ujuj) = (1.98 ± 0.08 (stat) ± 
0.32 (syst)) X 10^^. The corresponding product branching fractions with intermedi- 
ate resonances /o(1710), /2(1640), and /2(1910) are also determined, but with larger 
errors. 
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